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ABSTRACT 

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with 
a misaligned accretion disc to study BH spin precession and alignment jointly with 
BH mass Mbh and spin parameter a evolution, under the assumption that the disc 
is continually fed, in its outer region, by matter with angular momentum fixed on a 
given direction Jdiscout- We develop an iterative scheme based on the adiabatic ap- 
proximation to study the BH-disc coevolution: in this approach, the accretion disc 
transits through a sequence of quasi-steady warped states (Bardeen-Petterson effect) 
and interacts with the BH until the spin Jbh aligns with Jdisc.out- For a BH aligning 
with a co-rotating disc, the fractional increase in mass is typically less than a few 
percent, while the spin modulus can increase up to a few tens of percent. The align- 
ment timescale t a i oc a 5 / 7 M~ 32 / 35 is of ~ 10 5 — 10 6 yr for a maximally rotating BH 
accreting at the Eddington rate. BH-disc alignment from an initially counter-rotating 
disc tends to be more efficient compared to the specular co-rotating case due to the 
asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and 
opposite angular momentum when crossing the innermost stable orbit, so that the 
spin modulus decreases faster and so the relative inclination angle. 

Key words: accretion, accretion disc - black hole physics - galaxies: active, evolution 
- quasars: general 



1 INTRODUCTION 

Astrophysical black holes (BHs) are Kerr black holes fully 
characterized by their mass Mbh and spin Jbh, customarily 
expressed in terms of the dimensionless spin parameter a 
1), and unit vector Jbh: 



a Jbh- 



(1) 



The spin and mass of BHs residing in galaxy nuclei do not 
remain constant, close to their birth values, but change size- 
ably through cosmic time, in response to major accretion 
events. In current cosmological scenarios for the evolution 
of galaxies, repeated interactions among gas-rich halos play 
a key role not only in shaping galaxies, but also in trigger- 



ing quasar activity (White & Rees 1978 Di Matteo et al 



2005). Massive gaseous nuclear discs that form in the after- 



math of major galaxy mergers (Mihos & Hernquist 1996 



Mayer et al.]2007 l may provide enough fuel to feed, on sub 



parsec scales, the BH through a Keplerian accretion disc 
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1 Dotti et al.|2007 20091. If these episodes repeat recursively 
and/or at random phases (King & Pringle 20061 the BH 



spin Jbh is expected, initially, to be misaligned relative to 
the direction of the angular momentum of the disc Jdisc.out 
at its unperturbed, outer edge .Rout- In this configuration, 
the gas elements inside the disc undergo Lense-Thirring pre- 
cession (see, e.g. Wilkins 19721. In the fluid, the action of 
viscosity onto the differentially precessing disc ensures that 
the inner portion of the accretion disc aligns (or anti-aligns) 
its orbital angular momentum with the BH spin Jbh, out 
to a transition radius 7? W arp beyond which the disc remains 



aligned to the outer d isc, as first shown by|Bardeen fc Pet- 



terson 



(1975 1 (see also Armitage fc Natarajan||1999 



fc Papaloizou||2000[ |Fragile fc Anninos||2005| |Fragile et al 



Nelson 



20071. Warping of the inner disc at distance R from the BH 



is communicated through the fluid elements on a timescale 
fep(ii) related to the vertical shear viscosity of the accre- 
tion disc. Therefore, the inner regions of the disc align (or 
counter-align if the disc is counter-rotating) with the BH 
spin on the scale tfip(-Rwarp) when the viscous time for ver- 
tical propagation of disturbances equals the Lense-Thirring 
precession time. On a longer timescale, the joint evolution 
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of BH+disc system restores full axisymmetry, with the BH 
spin direction aligned relative to the total angular momen- 
tum of the composite system ( Rees|1978 Thorne et al.|1986 



King et al.|2005 I. The change in Jbh is a consequence of an- 



gular momentum conservation: since the BH acts on the disc 
with a torque that warps the disc, then an equal and oppo- 
site gravito-magnetic torque acts on the BH that modifies 
its direction only. 

BH spin alignment has been studied in two main con- 
texts. In the first, explored by |King et al.| ( [2005| and |Lodato| 
& Pringle ( 2006 1 , the focus is on an closed system where 
the accretion disc has a finite mass and radial extent. Here, 
the total angular momentum J tot = Jbh + Jdisc is well de- 
fined vector, and the BH eventually aligns it spin vector 
to the direction of Jtot. In the second, explored by |Scheuer 



fc Feiler| (|1996[ >, |Natarajan fc Armitage| ( |1999[ ) and |Martin" 



et al. (20071, the focus is on an open system, where the ac- 
cretion disc has infinite extension and it is continually fed at 
its outer edge by matter whose angular momentum has con- 
stant direction Jdisc, out- In this second case, the BH aligns 
its spin to the outer disc direction Jdisc, out on a timescale 
i a i that exceeds tBP by a few orders of magnitude ( jScheuer] 
fc Feiler|1996}|Martin et al.||2007 l. 



In this paper, we progress on the study of BH align- 
ment including the contemporary change in mass and spin 
modulus due to accretion of matter, neglected in previous 
works. During BH precession and alignment, matter flows 
inward and accretes carrying the energy and the specific 
angular momentum of the innermost stable circular orbit 
(ISO). This study thus provides estimates of the fractional 
increase of mass AMbh/A^bh,o an d spin Aa/ao during BH 
alignment (the subscript refers to initial conditions), to- 
gether with a sensible expression for the alignment time £ a i- 
In our context, we assume a continuous and coherent feed- 
ing of the accretion disc around the BH, at least for a time 
as long as the alignment timescale t a i. Thus, we consider 
an open system and we fix the orbital angular momentum 
direction Jdisc, out at the outer edge of the disc. 

In Section 2, we introduce key parameters and high- 
light our model assumptions. Section 3 surveys properties 
of steady state warped discs and key scales associated to 
the Bardeen-Petterson effect; disc models with constant and 
power-law viscosity profile are explored for completeness. In 
Section 4, we describe the equations for the BH mass and 
spin evolution, and introduce the adiabatic approximation 
to solve these equations. In the same section we also revisit 
the expression for the BH alignment time. Results are illus- 
trated in Section 5; there we explore also the tendency to 
alignment in initially counter-rotating warped discs. Section 
6 contains the discussion of the results and our conclusions. 



2 INITIAL ASSUMPTIONS AND MAIN 
PARAMETERS 

We consider a BH with spin Jbh, surrounded by a geomet- 



rically thin, standard Shakura-Sunyaev a-disc (e.g. Shakura 
|fc Syunyaev|1973| [Frank et al.|200"2| . The a-disc is initially 
misaligned relative to Jbh, i.e. the angular momentum unit 
vector of the disc at the outer edge is Jdisc, out 7^ Jbh; the 
relative inclination angle between the two unit vectors is 



a 


Ji/2 


0.18 


1.00 


0.15 


0.85 


0.09 


0.60 


0.05 


0.38 



Table 1. Table of the coefficients a and /„ 



Following |Pringle| ( |1992[ ) , we assume that the accretion 
disc has a high viscosity (a > H/R, where H is the disc 
vertical scale height) so that perturbations propagate diffu- 
sively. We introduce two viscosity parameters, v\ and v\ 
is the standard radial shear viscosity while v% is the vertical 
share viscosity associated to the diffusion of vertical warps 
through the disc, due to Lense-Thirring precession. For v\ 
we adopt the a prescription 



v\ = aHcs 



(2) 



where c s is the sound speed inside the accretion disc. It is 
still poorly understood which is the relation between the 
radial and the vertical viscosity; in particular, if v\ ~ V2 or 
v\ <tCvi- In order to simplify our discussion, we refer to the 
recent analysis of Lodato & Pringle (20071, and for V2 we 
take: 



^2 



f"2 

2a 2 



(3) 



where f V2 (given in Table 1) is a coefficient determined in 
numerical simulations that accounts for non-linear effects. 

The disc model is defined after specifying five free pa- 
rameters (subscript will be introduced to indicate initial 
values when mass and spin evolution is considered): 

(1) The BH mass, Mbh; we explore a mass range between 
1O 5 M < Mbh < 1O 7 M . For the BH mass we introduce 
the dimensionless parameter Me as Mbh = M D x 1O 6 M0. 

(2) The spin modulus, in terms of the dimensionless spin 
parameter a, which varies between < a ^ 0.95. We do 
not use the theoretical limit a = 1 because, if accretion 
is driven by magneto-rotational instabilities in a relativis- 
tic MHD disc, the final equilibrium spin due to continuous 
accretion is a ~ 0.95 (Gammie et al. 112004 1; 



(3) The relative inclination angle 9 out , between the spin ver- 
sor Jbh and the orbital angular momentum versor at the ex- 
ternal edge of the accretion disc, Jdisc, out- This angle varies 
isotropically from to n. In the following, however, we will 
confine this interval to (0, ~ 7r /6) in order to satisfy the used 
approximations . 

(4) The viscosity parameter a which is assumed to vary be 



tween 10 2 < a < 10 1 to bracket uncertainties ( King et al 



2007 1 . For our purposes we selected values of a according to 
Lodato & Pringle| |2007| , as in Table [I] In this study, a is 



considered as a constant inside the disc. 
(5) The accretion rate onto the BH, M, is expressed in terms 
of the Eddington ratio /Edd = L/L^dd and of the accretion 
efficiency 77 (where L^dd is the Eddington luminosity): M — 
/Edd-^Edd/^c 2 ). We consider values of /Edd in the interval 
10~ 4 < /Edd < 1 and compute r\ as a function of the BH 
spin modulus. 

If the disc, warped in its innermost parts, is described 
to first order by the Shakura-Sunyaev a-model, both v\ and 
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V2 follow a power-law. If viscosities satisfy relation Q and 
a is assumed to be constant, their exponent are equal: 



A V1 R" 



and 



A-vn Ft 







(4) 



Following standard Shakura-Sunyaev disc solutions for ex- 



ternal regions of an accretion disc (Frank et al. 2002), we 
have (3 = 3/4 and 



3/10 



5/4 -1 

cm ' s 



A V1 = 9.14 x lO^X 720 ( 

\ ^o.i 

A V2 = I I A V1 — 50 /„ 2 a ,i A V1 . 

In equation JSl, ao.i and 770.1 are the a coefficient and the 
BH radiative efficiency in unit of 0.1. f V2 is tabulated in 
Table 1 (Lodato & Pringle 2007). 



3 WARPED ACCRETION DISC 

3.1 The angular momentum content of discs: 
extended versus truncated discs 

The dynamics of a fluid element in a misaligned disc around 
a spinning BH is given by the combination of three different 
motions: the Keplerian rotation around the BH; the radial 
drift, due to radial shear viscosity, and finally the Lense- 
Thirring precession, due to the gravitomagnetic field H g 



generated by Jbh (see, e.g., Weinberg 1972] Thorne et al. 
|1986[ ). In response to Lense-Thirring induced precession, vis- 
cous stresses in the disc acts rapidly to produce in the vicin- 
ity of the BH an axisymmetric configuration whereby ad- 
jacent fluid elements rotates in the equatorial plane of the 
spinning BH. The disc thus warps and the warp disturbance 



propagates diffusely (Papaloizou & Pringle 1983 1 in the disc 



As the Bardeen-Petterson effect modifies the inclination 
of the orbital plane of consecutive infinitesimal rings, then 
the warped profile of the accretion disc can be described by 
the specific angular momentum density, L, expressed as 



L = LI = EQ. K R I 



(6) 



where l(R) is a unit vector indicating the local direction of 
the orbital angular momentum, L is the modulus, E is the 
surface density of the disc and Qk the local Keplerian an- 
gular velocity. The angle describing the tilted disc is defined 



0(7?) = cos~ 1 (l(i?)-J B H), 



(7) 



so that \(R) carries information of the warped structure of 
the accretion disc. The angular momentum of the accretion 
disc within radius R is given by 

Jdisc(-R) = f 2ivxL(x)dx (8) 
•> Riso 

where the integration domain extends from the innermost 
stable orbit -Riso out to R. In order to calculate the to- 
tal disc angular momentum we define an outermost radius, 
Rout - For an extended disc with i? ou t — > 00, the disc angular 
momentum Jdisc always dominates over Jbh. 

Real discs are likely to be truncated by their own self- 
gravity that becomes important at distances where the disc 
mass Mdisc^g) ~ (H/R) M B n (see, e.g., |Pringle|T98l] [Frank] 
et al.|2002 Lodato|2007 1. Outside the truncation radius, gas 



can be either turned into stars or expelled by winds from 



stars which do form (Levin 2007 King & Pringle 20071 



Thus, we are led to define a disc outer edge as the distance 
where the Toomre parameter for stability, Q = nc s / (nCE) 
(where k 2 — R(dQ 2 /dR) + 4f2 2 ), becomes less than unity, 
and the cooling timescale of the clumping gas is less than 
its dynamical timescale. When the Toomre parameter drops 
toward unity, the disc becomes unstable on a lenghtscal e 
A = c 2 /(GE) ( |Polyachenko et al.|l997||Levine et ah||2008| ); 
for a nearly Keplerian, Shakura-Sunyaev a-disc, this scale is 
much smaller than the disc radial dimension, and the cooling 
time of the associated perturbation is less or of the same 
order of its orbital period. Then, as long as the accretion disc 
can be described as a Shakura-Sunyaev discQ the external 
radius can be defined from the condition Q(R ou t) = 1, so 
that 

-22/45 



R out = 1.21 x W'afl i5 M { 



Rs, 



(9) 



-52/45 / /fidd 
V V0.1 

where Rs = 2GA/bh/c 2 is the Schwarschild radius. At the 
outer edge of the disc, l(R ou t) = Jdisc, out and 0(R O ut) = out - 
Definitions (JsJ and (JsJ) for L and Jdisc (R) hold for any 
disc profile. At first order, we can neglect details about the 
warped disc structure around R! W ar P assuming 1 « (0,0,1), 
and estimate the modulus of the orbital angular momentum 
within radius R, Jdisc (R), using Shakura-Sunyaev solutions 
for a flat disc. In this appr oximation, the surface density is 



Eflat ? 

and 
L(R) 



M I (37r^i ) (see, e.g. 



Pringle|l981 



Frank et 



al.|2002 l 



M 

37TJ/1 



VGM Bn R. 



(10) 



Using equations Q and ( 10 1, and espression Q for v\ in the 
case of f3 = 3/4, the modulus of the disc angular momentum 
within R reads: 



Jdisc (R) 



8 MVGMj 



BH^T/4 



21 



(11) 



If espression (111 is estimated at the outer radius |9|, the re- 
sulting dimensionless ratio between the disc and BH angular 
momenta is 

/ , \ -7/45 
-37/45 / /Edd 



Jdisc {R< 

Jbh 



out; _ „ 13/45 

= 7.3 a, 



0.1 M 6 



7/0.1 



(12) 



3.2 Timescales and warp radius 

The time-dependent evolution of the disc is described by the 
continuity equation 



r9 1 + 5s0*s*> = °. 



(13) 



where wr is the radial component of the velocity vector, 
and by the equation of conservation of angular momentum. 



1 This condition is fullfilled only for (M6/E,jd)/( a 0.l 170.1) ^ 4.3. 
If this condition is not satisfied the gas temperature drops below 
~ 10 4 K, in the external region of the disc where Q is still greater 
than unity. The change in the opacity likely modifies the structure 
of the outer disc, and we can not explicitly use (|9]l. In this paper 
we assume that the outer region is sufficiently extended to provide 
matter and angular momentum to the inner regions and use self- 
consistently the Shakura-Sunyaev model to describe the disc in 
regions where the gravitomagnetic interaction takes place. 
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In presence of a gravitomagnetic field, for a geometrically 
thin disc characterized by the two viscosities v\ and V2, the 
equation reads ( Pringle|1992 I: 



dt 



1 8 , Id 

RdR {RLVK)+ RM 



+ 



+ 



1 Br 91 



+ 



2G* Jbh x L 



(14) 



IP 



The last term is the Lense-Thirring precession term and the 
associated angular velocity is 

2G Jbh 
c2"~R3~' 



(15) 



The time-dependent equation ( 14 1 describes the radial drift 



of matter and the diffusion of warping disturbances across 
the high-viscosity disc. 

This equation introduces several key scales: 

(i) The viscous/accretion timescale for radial drift, re- 
lated to the angular momentum transport parallel to 
Jdisc,out, t llcc (R). It can be seen as the time it takes for a fluid 
element at R to accrete onto the BH (see, e.g., Pringle|1981 1 
Considering equation (141, the balance between the advec- 



tion term and the viscous term proportional to v\ (both on 
the right side of equation 14 1 leads to an estimate of the 
accretion time: 



R 2 / Vl . 



(16) 



According to equation (16 1, we can introduce the disc con- 
sumption timescale idisc, a concept useful when considering 
transient, truncated discs, as the accretion timescale at the 
outer radius: 

^disc ~ ^acc(-Rout ) 



1.71 x 10 6 a~i /45 M 6 



11/45 / /] 



-41/45 



(17) 



yr. 



(ii) The timescale for warp propagation, related to the 
radial diffusion of gravitomagnetic perturbations that trans- 
port the component of the disc angular momentum lying in 
the plane of the disc; this scale is inferred from equation ( 141 
considering the term proportional to V2, 

tnv(R) ~ — ~ (-) t«e(fl). (18) 



The physical interpretation of this timescale has been re- 



cently investigated by solving numerically equation ( 14 1 for 
a thin disc (Lodato & Pringle 2006): starting at t = with a 
flat disc misaligned relative to the fixed BH spin, t w 4bp (R) 
indicates the time it takes for the radial diffusion of the warp 
to reach radius -R; on longer timescale, the disc approaches 
a steady warped state. 

(iii) The characteristic extension of the warp -R warp , 
defined as the distance at which the Bardeen-Petterson 
timescale £bp(-R) equals the Lense-Thirring precession 
timescale O.^: 



Rw 



4GJbh 
V2C 2 



(19) 



For power-law viscosity model, equations |T]), (|5| and (191 

give 

-6/35 



24/35 
*0.1 



/,7 2 4/7 M 6 4/35 



no a 



a 4/7 R S . (20) 



The warp radius represents the dividing between the outer 



region for R ^> i?„ 



where the disc keeps its origi- 



nal inclination, given by Jdiscout, and the inner region for 
R <C -Rwarp, where the disc aligns (or anti-aligns) its orbital 
angular momentum with the BH spin, 1 || Jbh- The warp 
radius fixes also the magnitude of the relevant Bardeen- 
Petterson timescale, which reads 



iBp(^warp) — 33.5 a j /„ 2 i/. fi 



h 



72/35 f-12/7 ,,47/35 
18/35 



If we define the function 



5/7 

a yr. 



V(-R) 



r/1 

dR 



(21) 



(22) 



and Rbp the radius where the disc is maximally deformed 



* = V(-Rbp) = max(^) , 
we expect that: 

RbP = TlBpRwurp 



(23) 



(24) 



with ubp of order unity. -Rbp has two important proper- 
ties: first, if it is the radius where the disc is maximally 
warped, i.e. where the diffusive propagation of vertical per- 
turbations is more significant; second, it provides a reliable 
estimate of the distance from the BH where the gravito- 



magnetic interaction is stronger. From equation ( 14 1 this 
interaction is proportional to (L x Jbh)/^? 3 : this term van- 
ishes in the inner part of the disc (R <C -Rbp) since the 
Bardeen-Petterson effect aligns L with Jbh, and also in the 
outer regions (R 2> -Rbp), due to the rapid decline with R. 
Accordingly, the region near Rbp (or equivalently -R W ar P ) is 
the only one significantly misaligned with Jbh. 



3.3 Analytical solutions 

In this Section we summarise the properties of the steady 
warped disc structure used to compute the joint evolution 
of the disc and the BH. 

Following previous studies we assume that the viscosity 
profiles are power-laws with exponent f3, as in equation Q, 
and explore two possible cases. In the first, we formally ex- 
tend the Shakura-Sunyaev solution everywhere in the disc, 
i.e. v\ oc R? with (3 = 3/4, and V2 given by equation ^ 
( jMartin et al. 12007] ). In the second case, we assume the vis- 
cosities to remain approximately constant everywhere in the 
disc ( |Scheuer fc Feiler||1996[ ). In order to compare the two 
models ( jMartin et al.||2007| cfr.), we impose the continuity 
of the viscosities at Rbp where the gravitomagnetic torque 
is most important. 

Before solving equations ( |13[ ) and ( |14[ | we introduce two ap- 
propriate reference frames. The first is the inertial reference 
frame Oxyz referred to the outer disc; we can always rotate 
it so that its 2 axis is parallel to the direction of Jdisc,out- 
The second reference frame is the not-inertial frame O'x'y'z' 
referred to the BH spin, which is always centered to the BH 
and whose z axis is parallel to the black hole time varying 
spin Jbh ■ If we use the adiabatic approximation, then frame 
O'x'y'z' can be approximated, with time t, as a sequence of 
frames, one for every quasi-stationary state of the system. 
The shape of the warped accretion disc is studied in the 
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Figure 1. ip as function of i?/-R W arp, for Wq (solid line) and 
Wp L (dashed line) for different inclination angles: black lines 
correspond to 9 ou t = 7r/3, blue lines to n/30, red lines to 
7r/300. The vertical black dotted line represents Bep/Kwarp = 
0.42, the Bardeen-Petterson radius for constant viscosity pro- 
files. The parameters set for the BH and the disc is given by 
M BH ,o = 1O 6 M , a = 0.5, / Edd = 0.1, a = 0.09 



Figure 2. Xp—o{R) (solid line) and X/3=3/4 (dashed lines) as 
functions of i?/i? W arp for different inclination angles. Black 
lines correspond to 8 ou t = 7t/3, blue ones to 7r/30, red ones to 
7r/300; the vertical black dotted line represents i?Bp/^warp = 
0.42, the Bardeen-Petterson radius for constant viscosity pro- 
files. The parameters set for the BH and the disc is given by 
M BH ,o = 1O 6 M , a = 0.5, / Edd = 0.1, a = 0.09 



O'x'y'z' frames and the cartesian components of any vector 
v are there indicated as v' x , v' y , v' z ; Oxyz is the natural frame 
to study the temporal evolution of the BH spin and here the 
Cartesian components of the previous vector are denoted as 



For a stationary state, continuity equation ( 13 1 can be easily 



intergrated introducing the accretion rate M as constant of 
integration: 



RT,v R 



M 



while the projection of equation ( 14 1 along 1 reads: 

= 



3 dL MVGM^\ , 1_ _ 
2 dR 4 n ^R J 2 



dl 



dR 



(25) 



(26) 



equation ( 14 1 in the small deformation approximation have 
been found by|Martin et al.| pOOT]): 



B 



x K x 



R 



Rw 



/2(l+/3) 



^2(1 



(1 + 0) 



R 



Rv 



(28) 



where B is a complex constant of integration, depend- 
ing on the boundary condition at the external edge, the 
subscript "PL" is a reminder of the power-law viscosities 



and K ± 



/(2(l+/3)) 



is the modified Bessel function of order 



1/(2(1 + /3)). In the particular case where we consider con- 
stant viscosities, i.e. /3 = 0, the solution can be written as 



In the small deformation approximation JScheuer & Feiler ^ G ^4 ex P I v2 (1 i) 



1996 ) the warp is gradual and we can neglect the non-linear 
term, proportional to \d\/dR\ 2 . Using the boundary condi- 



R 



Rw 



(29) 



tion E(7?iso) = 0, the integral of (261 is 
M 



L{R) = 



37T^i 



VGMbhR 1 



Riso \ 
R 



(27) 



This means that, in this approximation scheme, the modu- 
lus of the angular momentum density for a warped accretion 
disc far from the horizon is the same as for a flat disc, equa- 



tion ( 10 1 



Following Scheuer & Feiler ( 1996 1 we study the disc 



profile of the steady disc introducing the complex variable 
W' = l' x + il'y and considering the case O ut < tt/2. Using 
power-law viscosities according to Q, analytic solutions of 



where A is a complex constant of integration and the sub- 
script "C" is a reminder of the constant viscosity model 
I Scheuer fc Feiler|1996 1. In this latter case, ubp is calculated 
self-consistently, using the definition ( |24[ ) and the prescrip- 
tion for constant viscosity evaluation at i?p,p; we find that, 
for every possible parameters set, ubp ~ 0.42. We notice 
that W' (and so also ip) depends on the radius 7? through 
the dimentionless ratio R/R ws . rp ( Martin et al.|2007 1 . 

In Figure [l] we plot the modulus of the gradient of i, 
tp{R), which is a local measure of the deformation degree 
of the disc, for a particular set of parameters and for three 
different angles, (9 out = 7r/3, 7r/30, 7t/300. The shape of if) is 
similar for the two different disc profiles and for all the an- 
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expect the disc to be more deformed. At radii smaller than 
^warp and far from -R W ar P the disc is almost flat (note that 
the graph is logarithmic in both axes). For the constant vis- 
cosity (power-law) profile the peak is at _Rbp ~ 0.42i? warp 
(Rbp ~ 0.38-Rwarp). In Figure 1 we also see that a constant 
viscosity disc is less warped (since the maximum deforma- 
tion is the smaller) than the power law viscosity disc. The 
ratio between the maximum deformations in the power law 
vs constant viscosity is roughly a factor 2, and it does not 
depend on the inclination angle (except for a scale factor, 
approximately equal to the ratio between the corresponding 
angles) . 

3.4 Validity of the approximation 

We calculated the warped disc profile under the small defor- 
mation approximation. We neglected second order terms in 



equation (26 1 and found an analytic solution for L; in order 



to verify the consistence of this approximation, we define X/3 
as the ratio between the neglected term and the first term 



into the round brackets of (261, assuming to have a Keple- 



rian disc with power-law viscosity profile with exponent j3, 
like in equation Q. Considering equations (271 for L, we 
have dL/dR « (1/2 + /3)(L/R) and then \P reads 



xp(R) 



_rfl_ 

Jr 



R z 



(30) 



3 (/? + 1) pi 

Once we know the explicit solutions, the consistence of this 
approximation can be tested a posteriori calculating \P- the 
approximation is well satisfied if Xf3 <S !• From equation 
( 30 1, X0 can b e expressed also as a function of R/ Rw&rp and 

2 Vn „9 / R \ 2 ,0. ( R 



- — R 2 

X0 ~ 3 09 + 1) vi p 



r 



(31) 



Figure [2] shows the function X/3=o for constant viscosity pro- 
files (dashed lines) and the function X/3=3/4 for power-law 
viscosity profiles (solid lines), for the same parameters as in 
Figure [l] 

The function Xs exhibits a maximum, (x/3)max, around 
fl warp . Far from 7? W ar P the accuracy of the approximation 
increases, albeit slowly. The function xp is most sensitive to 
the inclination angle, as expected (notice that Figure 2 uses 
logarithmic axes). 

In Figure [3] we test the validity of the small deforma- 
tion approximation plotting, in the BH mass versus 6 ou t 
plane, the color coded values of (x/3=o)max, for different val- 
ues of the viscosity parameter a (Table [T| , using the con- 
stant viscosity profile model (we fix /Edd = 0.1 and a = 0.9). 
White zones represent the regions where (x/3=o)max > 1, 
i.e. where the small deformation approximation becomes in- 
valid. (x/3=o)max shows mainly a strong dependence on in- 
clination angle # ou t, but also a weaker dependence on the 
BH mass which reveals that the small deformation approxi- 
mation is less accurate for Mbh > 10 6 Mq and increasing 
BH mass. Comparing different a values, the approxima- 
tion is better satisfied for large viscosities parameters (i.e. 
a — 0.18). We repeated the analysis for the power-law vis- 
cosity model that shows no significant differences in the pa- 
rameters dependence. 

In Figure[4j using the same colours conventions, we explored 
(X/3)max in the ou t versus a (left panels), and /Edd versus 



#out (right panels) planes, once we have fixed the viscosity 
parameter (a = 0.09), the BH mass (M B h = 10 6 M Q ), and 
/Edd = 0.1 for the left panels and a = 0.9 for the right pan- 
els. For both constant (f3 = 0) and power-law (J3 = 3/4) 
models the relative inclination angle is again the leading pa- 
rameter gauging the goodness of the fit as the approximation 
depends very weakly on a and /Edd- 



4 BLACK HOLE EVOLUTION 
4.1 Basic equations 

In this section we explore the equations for the BH evolution. 
The BH is accreting and its mass increases, from an initial 
value Mbh.o, according to 



dMn 



tit 



M 



E(Rd 



(32) 



where E(Riso) is the energy per unit mass of a test particle 
at the innermost stable orbit. E(Rjso) / '<? = 1 — v( a ) is 
related to the efficiency 77(a) that depends only on the spin 



parameter ( |Bardeen 1970 Bardeen et al. 19721. Equation 
( 32 1 introduces a natural timescale for BH mass growth, 
known as Salpeter time ts: 



ts = 4.5 x 10 f 



/Edd(l - rj) 



yr. 



(33) 



As argued by Rees ( 1978 1 and shown by Thorne et al 



( 1986 1, there is a coupling between the BH spin and the an- 



gular momentum of the disc. Even though the disc is much 
less massive that the BH, the moving fluid elements perturb 
the Kerr metric and interact with the BH spin, causing spin 
precession, and if viscous dissipation is present, alignment. 
For an infinitesimal ring of inviscid matter with total angu- 
lar momentum Jring, the BH spin precesses, following the 
equation 



dJ BH 2G Jrin 



X J 



BH, 



dt C 2 R 3 

with a precession frequency 



/-jprcccssion 

i ^BH — 1 



J ring 
^BH 



(34) 



(35) 



Equation ( 34 1 can be extended to the case of an accretion 
disc to yield: 



dJ E 



dt 



MA(R 1S o)l{Rv, 



4nG 



L(R) x J BH 
R? 



(36) 



dR. 



The first contribution is due to accretion of matter at J?iso 
where A(_Riso) indicates the orbital angular momentum per 
unit mass carried by matter at ISO; the Bardeen-Petterson 
effect ensures that the direction of I(i?iso) is parallel or 
anti-parallel to Jbh, so that the accretion modifies only 
the spin modulus. As shown by |Bardeeri] ( |1970| ), a varia- 
tion of mass AMbh = v6Mbh,o is necessary to pass from a 
Schwarzschild BH (a = 0) to an extreme Kerr BH (o = 1), 
while spin flip of n, due only to accretion on an initially 
extreme Kerr BH, needs AMbh = 3A/bh,o- So, the spin 
accretion timescale for the spin modulus is of the same or- 
der of the mass accretion timescale ts- The second term 



in equation ( 36 1 describes the gravitomagnetic interaction 
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Figure 3. Color coded plot of (xfl)max in the 8 ou t versus Mbh plane, for four different a parameters: a = 0.05 top left panel, a = 0.09 
top right panel, a = 0.15 bottom left panel, a = 0.18 bottom right panel. The disc has constant viscosity profiles, i.e. (3 = 0. The 
accretion rate is /Bdd = 0.1 and the spin parameter is a = 0.9 



between the rotating viscous disc and the BH spin vector. 
This term modifies only the spin direction of the BH in or- 
der to conserve the total angular momentum of the system. 
Under the working hypothesis that the disc is continually 
fed by matter carrying the same angular momentum (see 
Section [6] for a critical discussion) , the BH aligns its spin 
Jbh in the direction of Jdisc.out- Alignment implies that 
0out(t) = cos _1 (Jbh(*)- Jdisc.out) goes toOwithtime. Figure 
[5]shows the function 7 defined as the modulus of the integral 
kernel of equation ( 36 1 



7(7?) = 



47tGL(7?)Jbh bw.[6(R)] 



c 2 R 2 
as a function of 7?/7? warp , 



(37) 



7r/3, 7r/30, 7r/300, where 6(R) is computed along the pro- 
file of the steady warped disc of equation i\28\ and ( 29 1 . The 
function 7, similarly to ip (defined in eq. |22|), peaks near 
7?warp. Contrary to ip, power law viscosity profiles have lower 
peaks, compared with constant viscosity profiles. This figure 
indicates also that the BH-disc gravitomagnetic interaction 
is spread over a relatively small region of the disc around 
the warp radius; the characteristic spreading length, which 
is slightly larger for constant viscosity profiles, is usually of 
a few warp radii. 
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Figure 4. Left panels: color coded plot of (x^)max in the a versus 8 ou t plane, for a = 0.09, M BH = 1O 5 M and / Edd = 0.1. Right 
panels: color coded plot of (X/3)max in the /Edd versus 9 ou t plane, for a = 0.09, Mbh = 10 5 Mq and a = 0.9. Top panels refer to constant 
viscosity profiles, bottom panels to power-law viscoisty profiles. 



4.2 Alignment time 



and the integral on the right hand side as 



In this Paragraph we want to give simple estimations for 
the alignment and the precession timescales, starting from 



equation ( 36 1 



Assuming BH mass and spin modulus variations due to 
accretion to be small compared with gravitomagnetic effects 
during the alignment, we neglect the term proportional to 



A(-Riso) in (36 1; if BH spin aligns and precess, left hand 
side of ( 36 1 can be estimated introducing a characteristic 
gravitomagnetic timescale r gm as 



dt 



A Jbh I Jbh sin # O ut,0 



4ttG 



L(g) x Jbh 



dR 



4-kG L{R warp)t/l3H sin o 
J 1-warD 



since the bulk of the gravitomagnetic interaction occurs 
around Rv, 



Equating these two expressions and using 
equation (271 for the specific angular momentum density 



modulus, we obtain 



3c^i(i? warp ) R 



GM 



(■warp 



(38) 



Using equation ( 19 1 
My/ GMbh /ui(R warp / 



and (111 which imply 

(21/8) Jdisc(-Rwarp) -Rwarp , the 



gravitomagnetic scale T gm can be written in terms of the 
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Figure 5. In this figure we draw modulus of gravitomagnctic 
interaction term, /, as a function of radius normalized to warp 
radius, i?/i? W arp- Disc profiles are obtained by BH with Mbh = 
1O 6 M0 and a = 0.5, and accretion disc with /Edd = 0-1 and a = 
0.09. Solid (dashed) lines refer to constant (power-law) viscosity 
profiles; black lines to f? out = tt/3, blue lines to 8 ut = tt/30, red 
lines to 6 ou t = 7r/300. 



Bardeen-Petterson warp timescale (eq. 18 
iV2 



Jbh 



-tep(^warp) 



7 Jdlsc(R warp } 

and also in term of the accretion timescale (eq. [l6] 

7 1*2 Jdisc(Rwarp)- 



iacc (^Vi 



(39) 



(40) 



where Jdisc(-Rwarp) is the disc angular momentum modulus 



within the warp radius, estimated by (111. Finally, consid- 
ering equation (111 and (161 for i aC c together with the ex- 



pression for the spin modulus and Schwarzschild radius, r g 



of expression ( 40 1 can be rearranged as 



3 v\ Mbh 

- a -. — 

2 u 2 M 



Rs 



R™ 



(41) 



Since the disc carries very little angular momentum at 
the warp radius, from equation (39 1 r gm ^> t B p, always. The 



gravitomagnetic BH-disc interaction causes BH spin preces- 
sion and alignment at the same time, and then introduces 
two scales related with r gm , the precession and the alignment 
timescales, t plcc and i a i respectively. We separate their rel- 



ative importance following Martin et al. (20071 results, and 
define the parameter fi, so that 



iai 



cos /j, 



sin ji 



(42) 



The exact value of /i depends on the viscosity profile, and can 
be estimated either analytically ( Martin et al.|2007| ), or nu- 
merically as in this paper. Initially, we assume alignment and 
precession to have the same timescale, cos/i = sin/i = V2/2 



according to Scheuer & Feiler (19961. Substituting expres- 
sions ^ for the viscosities, (20 1 for the warp radius, and 
(|38[) for r gm in (42 1, the alignment time reads 



1.13 x 10 q 



5^58/35 j-5/7 A/f -2/35 ( /Edd 



M,7 



VO.l 



-32/35 



5/7 

a yr. 



(43) 



The timescale i a i increases with o, indicating that a rapidly 
rotating Kerr BH offers some resistance before changing its 
direction. Interestingly, the alignment timescale does not de- 
pend on the initial inclination f? out ,o since a more inclined 
configuration implies more pronunced disc deformations and 
stronger mutual gravitomagnetic interactions (as also shown 
in Figure [I] and [5]) . i a i has a weak dependence on the BH 
mass and scales nearly as M _1 : a higher accretion rate im- 
plies a higher angular momentum density L(i?) and thus 
a stronger gravitomagnetic coupling. We notice also that, 
apart from numerical factors of order unity, this timescale 
is consistent with the alignment scales found by by|Scheuer 
fc Feiler| ([19961); |Natarajan fc Pringle] (|1998[); |Natarajan"& 



Armitagc 



< |1999| ); |Martin et aL] ( [2007 1 



4.3 The adiabatic approximation 



In Section |3.2| and |4.1| we described the equations govern- 
ing the evolution of a warped accretion disc around a fixed 
BH, and the evolution of an accreting Kerr BH in gravit- 
omagnetic interaction with its accretion disc. The BH and 
the accretion disc evolve contemporary and their evolution 
is coupled, so that we can solve simultaneously equations 



( 13 1 and ( 14 1 for a Keplerian disc and ( 32 I and ( 36 1 for the 



accreting and precessing BH. 

In this paper, we solve these coupled equations using 
the adiabatic approximation that separates the rapid tem- 
poral evolution of the warped disc from the longer temporal 
evolution of the BH. Equations are integrated starting from 
given initial conditions: at t = the BH spin Jbh is in- 
clined with respect to Jdisc,out by an angle O ut,o> an( i the 
warped disc profile is described by a quasy-stationary profile 
L(R,t = 0); A/bh,o and Jbh.o are the initial BH mass and 
spin. 

In order to justify this approximation scheme, we survey 
the BH and disc timescales, as functions of Mbh and /Edd, 
for two selected values of the viscosity and spin parameter: 
a — 0.15 and a — 0.9. In Figure Q and in Figure Q, we 
draw in the MBH-/Edd plane lines of constant tBp(R w )/t il i 
and ts/tsd ratios. The comparison between the different 
timescales lead to the following hierarchy of timescales: 



t B p(R 



< tai < iS- 



(44) 



Then, in the adiabatic approximation, the disc transits 
through a sequence of warped states over the shortest 
timescale iBp(-Rwarp) while, on the longer timescale i a i, the 
BH aligns its spin to Jdisc.out, and modifies a little its spin 
modulus and mass due to accretion. Considering one of these 
disc quasi-steady states, initially at time i, after a time gap 
St ~ iBp(i?war P ) the BH mass and spin Jbh are updated 
according to 



M BH (t + iBp(i?war P )) = Mbh (i) + 5M B n 
Jbh (i + iBP (-Rwarp ) ) = Jbh (i) + 5 Jbh 



(45) 
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Figure 6. In the Mbh versus /Edd plane, we draw lines of 
constant tBp(^?warp)/i a l ratio for a BH with a = 0.9 and an 
accretion disc with a = 0.09. 



Figure 7. In the Mbh versus /Edd plane, we draw lines of 
constant ts/i a l ratio for a BH with a = 0.9 and an accretion 
disc with a = 0.09. 



and these variations produce a new quasi-stationary warped 

State at t + tBp(-Rwarp), L(i?, t + iBp(-Rwarp))- 

For the BH mass variation 8 Mbh, we integrate equation ( 32 \ 
from t to t + tpp(RBp)' 



5M B 



Ki E(Rko) tBp(RBp) 
Cr 



(46) 



where 7trso is the last innermost stable orbit associated with 
the current value of a(t). 

For the spin variation, we need to integrate equation ( 36 1 
that includes the two different and coupled contributions 
due to accretion and gravitomagnetic interaction; if <5Mbh 
and <5Jbh are small on the timescale tBp(RBp), to first order 
the two contributions decouple and they can be integrated 
separately: 



(SJb 



(Mi 



MA(Riso) t BP (RBp) 



=— tBp(RBp) 

C 

L(R,t) x J B h(*) 

disc 



(47) 



(48) 



R 2 



dR 



where (5JBH) acc is due to accretion and changes only the 
spin modulus while (<5JBH) gm is due to gravitomagnetic in- 
teraction and changes only the spin direction. After the in- 
terval iBp(flW), the angular momentum of ( ]45| are updated 
according to this rule 



(49) 



Jbh(£ + tBp(-Rwarp)) = ( Jbh(*) + (<5Jbh) 

„ JbhW + (<5J BH ) acc 



Jsn(t) 

This procedure can be repeated iteratively on a timescale 
ts\ to study the coupled evolution of ~L(R,t), Jbh and Mbh 
during the alignment process. 



5 SPIN ALIGNMENT 
5.1 Set up 

In this Section we study the coupled evolution of the BH 
and warped accretion disc using the approximation scheme 
described in the previous Section, in order to infer the 
evolution of Mbh and Jbh as a function of time, in response 
to the gravitomagnetic interaction and matter accretion. 

At t = the outer disc, extending up to a radius Rout, 
defines the fixed reference frame Oxyz. In this frame the 
external edge of the disc lies in the x, y plane and the orbital 
angular momentum at R ou t is 



(Lx (Rout ), L y (Rout), L z (Rout)) — L(R ou t)(0, 0, 1) 



(50) 



while the BH spin is initially inclined of #0,^0 with respect 
to the z axis: 



(Jbh, x, Jbh, v , Jbh, z) = JBn(sin (9 ou t,o, 0, cos # ou t,o). 



(51) 



If at t ^ we know the components of Jbh in the fixed ref- 
erence frame Oxyz there is always a rotated reference frame 
Ox'y'z' where Jbh is along the new z' axis (see also the 
discussion about reference frames of Section |3.3[ ) . The two 
reference frames are related by a rotation 1Z, which depends 
only on the components Jbh, a;, Jp.h, v , Jbh,z of Jbh(£) in the 
fixed reference frame. If TZij is the matrix associated with 
this rotation, we can easily find the components of Jbh(£) 
and li(R ou t,t) in the rotated frame: 



JknAt) =1l, 3 JBH,j(t) J BH (*) = (0,0, J BH (*)) 
Li(Rout,t) = TZij Lj (Rout, t)- 



(52) 



As shown by Scheuer & Feiler ( 1996 1 for the constant vis- 



cosity profile and Martin et al. (20071 for the power-law 



viscosity profile, in this special rotated frame of reference 
it is possible to calculate analytically the expression of the 
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Figure 8. Results for precession and alignment processes. Black lines refer to our result while red lines refer to results published by |Martin] 
ct al. ( 2007 ) ; solid lines (dashed lines) refer to constant (power-law) viscosity profile. Top left panel represents temporal evolution of relative 
inclination angle 9 ou t while top right shows evolution of Jbh,x/-'bh agains JBH,y/^BH> both for an initial BH with Mbh.o = 10 6 Mq, 
ao = 0.5 and an accretion disc with f^dd = 0.1 and a = 0.09, with O ut,o = 7r/6. Blue dashed line represents the evolution of the spin 
components for a pure precession motion around Jdisc.out II in bottom left (right) panel we represent evolution of Jbh x /Jbh agains 
>/BH,y/^BH for an initial relative inclination angle O ut,o = f/30 (# O ut,0 = t/3)j for an initial BH with Mbh = 1O 6 M0, ao = 0.5 and 
an accretion disc with /sdd = 0-1 an d Ct = 0.09. 



gravitomagnetic torque, using equation ( 48 1 



. 47rG Jbh (t) 



Jdis 



L(R,t) W'(R,t) 
R? 



dR) t B p(R 



(53) 



where L(R,t) is given by (27 1 and W'(R,t) by (29 1 or (281 



The analytic expressions of the gravitomagnetic torques for 
the two viscosity profiles are reported in the Appendix. From 
the torques it is possible to find the values of the spin varia- 
tions (<5JBH,ic)gm, (SJbh ,j/)gm and (5JBH,z)gm in this rotated 
reference frame. Finally, in order to know their expressions 
in our fixed reference frame Oxyz, we have to rotate them 
back, using the inverse rotation 1Z~ 

(<5JBH,i)gm = {R-~ 1 )ij($JBH,j)gm (54) 



Once we know the spin variations due to gravitomagnetic 
coupling, the modulus variation can be calculated from 



equation (471 and the global spin variation from equation 
(491. 



5.2 Results 

We computed, within the adiabatic approximation, the joint 
evolution of the BH mass and spin during the process of 
alignment under the assumption that matter is corotating 
with the BH. We iterated equations (1461) and (|49l), from 



the initial conditions (50 1 and (511, recording the updated 



values of Mbh, a and of the relative inclination angle 9 out 
every snapshot of time St ~ fBp(-Rwarp)- We initially choose 
a spinning BH with Mbh = 10 6 Mq and a = 0.5, and an ac- 
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Figure 9. Coupled evolution of the relative inclination angle 9 ou t, BH mass Mbh and spin parameter a. Solid (dashed) lines refer to 
constant (power-law) viscosity profile. Black lines to /Edd = 1; red lines to /Edd = 0.1, blue lines to /Edd = 0.01. Dotted horizontal lines 
which appear in top panels represent angles 9 ou t /0 O ut,o = 10 _1 , 10~ 2 , 10 — 3 . Initial configuration: Mbh.o = 10 6 Mq, ao = 0.5, /Edd = 0.1 
and a = 0.09, with O ut,o = t/6. 



cretion disc with M — 0.1 MEdd, a = 0.09; both power-law 
and constant viscosity profiles are considered. Three initial 
relative inclination angles of # O ut,0 = tt/3, 7r/6, 7r/30 have 
been tested. Figure [8] shows as a function of time the in- 
clination angle 8 ou t(t) and the two components of the BH 
spin unit vector in the plane Oxy; red lines refers to the an- 



alytic solutions given by |Martin et al. (20071. As shown in 
top-left panel of Figure 18] the relative inclination angle # out 
decreases exponentially with time on the scale tai and the 
decrease is more rapid for the constant viscosity disc (solid 
line) . 

The BH spin aligns with the external disc and precesses, 
as illustrated in the the top-right panel, where we compare 
also the actual evolution of the spin versor with a pure pre- 
cessional motion (blue dashed line). Our results are only 
qualitatively consistent with Martin's results; in our calcu- 
lations the alignment process appears to be less efficient 
and the spin precession more pronounced. The difference 



between our results and Martin's analytical solutions arises 
from three facts: (i) we included mass and spin modulus 
evolution; (ii) Martin et al. neglected to carry out the rota- 
tion connecting the BH reference frame O'x'y'z' to the disc 
frame Oxyz; (iii) for constant viscosity profile, we evaluate 
fi and V2 from equation |5| at the Bardeen-Petterson ra- 
dius, 7?bp ~ 0.4_R warp , while Martin et al. evaluate them at 
the warp radius, -R W ar P . For an initially not very inclined BH 
spin, the difference tends to disappear, because the rotation 
matrix nears the identity matrix. However, for large # O ut,0 
the discrepancy becomes more important (see, e.g., bottom 
panels of Figure [sj . 

Figure [9] shows the evolution of # ut and a as functions 
of time and of the increasing BH mass, for an initial BH with 
Mbh,o = 10 6 Mq, # O ut,0 = t/6 and spin parameter ao = 0.5, 
and for /Edd = 1, 0.1, and 0.01. Both constant and power-law 
viscosity profiles are explored, always with viscosity param- 
eter a = 0.09. Alignment is a process that shows a strong 
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IO" 1 


2.18 


3.9 


0.54 


2.30 


PL 


0.1 


10~ 2 


4.39 


7.9 


1.08 


4.57 






10~ 3 


6.64 


11.9 


1.64 


6.84 






io- 1 


18.1 


3.9 


0.45 


1.91 




0.01 


10~ 2 


36.3 


7.9 


0.89 


3.79 






IO" 3 


54.7 


11.9 


1.35 


5.67 



Table 2. Summary of our parameters and results for the co-rotating case; we consider viscosity coefficient a = 0.09 and initial inclination 
angle 9 out ,o = t/6, both for constant (C) and power-law (PL) viscosity profiles (VP). The initial BH has AfgH = 10 6 Mq and ao = 0.5. 
Accretion rate /sdd varies over three orders of magnitude and we record times needed to decrease the relative inclination angle of a 
factor 10, 100 or 1000, comparing it with estimated alignment timescale, equation (|43|; we also report mass and spin relative variations. 



dependence on the accretion rate: for the constant (power- 
law) viscosity model the time necessary to reduce the rela- 
tive inclination angle by a factor 100 varies from 3.0 x 10 5 yr 
(5.3 x 10 5 yr) for / Edd = 1 to 1.86 X 10 7 yr (3.63 x 10 7 yr) 
for /Edd = 0.01. During this alignment time, the BH has in- 
creased its mass by a small fraction, between 0.74% (1.30%) 
for / Edd = 1 and 0.46% (0.89%) for / Edd = 0.01. The spin 
parameter a increases due to accretion, but only by a small 
amount, between 3.13% (5.47%) for / Edd = 1 to 1.96% 
(3.79%) for / Edd = 0.01. 

In Table [2] we summarize the results of Figure [9] There, 
we compare also the time At necessary to decrease the ini- 
tial inclination angle by a given amount, with t a i estimated 



from equation (43 1: the values of the time At are consis- 



tent with the interpretation of i a i as e-folding time. For con- 
stant viscosity profiles, a closer match of t a i with the nu- 
merical outcomes requires (cos fx)c ~ 0.78 instead of y2/2, 
and (cos /x)pl ~ 0.41 for power-law viscosity profiles (sub- 
scripts C and PL simply remind that for different viscosity 
prescriptions we found different cos/x values). Then, the ra- 
tios between the precession and the alignment timescales 
are (t prC c/tai)c = 0.81 and (t prec /tai)p L = 2.2. These result s 
are still qualitatively consistent with Martin et al. (20071, 



who have shown that both t a j and t prec /t a i increase with 
the exponent /3 of the viscosity profile; small quantitative 
differences are due to different assumptions and different 
calculations methods. Finally, we notice that the scaling of 



At with /Edd is in good agreement with estimation ( 43 1 



5.3 Exploring the parameters space 

Here we explore more systematically how the fractional in- 
creases of Mbh and a, and the alignment time vary with ini- 
tial mass Mbh.o, spin ao, / E dd and a, for both constant and 
power-law viscosity profiles, fixing 8 ou t,o = 7r/6. The evolu- 
tion is followed until 9 out has decreased by a factor 100; we 
define as Atg _g / 100 the corresponding "alignment" time, 
computed self-consistently. We also infer from the numerical 
model the relative growths of BH mass AMbh/A/bh,o and 
spin parameter Aa/ao during At Bo ^g / 100 . 

Figure [lO] and Figure [Tl] show the weak dependence 



of the alignment time At, 



e ^e /ioo 



on the initial BH mass 



Mbh.o , and of the relative mass and spin parameter in- 
creases, for eight different sets of the other parameters. Com- 



paring numerical scaling factors for Me in Atg ^ 



e /ioo 



with 



that of expression (43 1, we notice again a good agreement, 



in particular for /Edd not too close to the Eddington limit 
and M BH ,o < 10 6 M Q . 

By contrast, the alignment process is more sensitive on 
/Edd, ao and a. Color-coded maps of At8 _>e /ioo (Figure 
[12] ), of AM/M BH ,o (Figure [13} and Aa/a (Figure [H]) are 
constructed in the ao versus /Edd plane, varying the coeffi- 
cient a and the viscosity law inside the accretion disc. 
In Figure [12] we infer the interval of the alignment time 

—>e /ioo (as inferred from the numerical model) of in- 
terest for the study of BH evolution. The alignment time 
can vary by many orders of magnitude from ~ 10 5 yr to 
~ 10 10 yr, and it reveals strong dependencies both on the 
accretion rate and on the initial spin parameter. In addi- 
tion, smaller viscosities (a = 0.09) gives shorter timescales 
compared to higher viscosities (a — 0.18). A simple com- 
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Figure 10. Alignment time (defined as the time needed for 
the relative inclination angle to reduce of two orders of mag- 
nitude, going from n/6 to 7r/600 ), as a function of the initial 
BH mass, Mbh o> f° r constant viscosity profile. The black lines 
refer to /Edd = 1 an d the red ones to /Edd = 0.001; the solid 
lines are for initial spin parameter ao = 0.9 while the dashed 
ones ao = 0.1; finally, the thin lines represent a = 0.09 and 
the thick ones a = 0.18. 



Figure 11. Mass and spin relative increase during the align- 
ment time (defined as the time needed for the relative in- 
clination angle to reduce of two orders of magnitude, go- 
ing from 7r/6 to 7r/600 ) as a function of the initial BH 
mass, for constant viscosity profile. Black (red) lines refer to 
/Edd = 1 (/Edd = 0.001). Solid (dashed) lines are for ao = 0.9 
(ao = 0.1); finally, the thin (thick) lines represent a = 0.09 
(a = 0.18). 



parison between alignment times Atg _g / 100 for different 
initial spin parameters, but identical /Edd, reveals that the 



scaling factors for a and 770.1(a) in equation (431 are in good 
agreement with numerical results. 

Figure [13] shows that the relative amount of mass accreted 
during the alignment process is small, compared to the ini- 
tial BH mass. It varies between ~ 10 -3 and ~ 10 -2 for the 
constant viscosity profile, and between ~ 2.5 x 10~ 3 and 
~ 3 x 10~ 2 for the power-law viscosity profile. Even if the 
accretion rate varies over four orders of magnitude, there are 
no comparable variations for the relative BH mass growth, 
in fact a larger /Edd means a larger accretion rate, but it 
also reduces the alignment time. The relative increase of the 
spin parameter a is shown in Figure [l4"| The evolution of a 
is the combination of different, and sometime opposite, ten- 
dencies: a highly spinning black hole requires a longer time 
to align, but the particles at its innermost stable orbit carry 
on the BH a smaller angular momentum. The spin modu- 
lus increases significantly during the alignment for initially 
slowly rotating BHs and high accretion rates, typically with 
5 x 10 -3 < Aa/ao < 8 x 10~ 2 for a constant viscosity profile 
and 10" 2 < Aa/ao < 2 x KT 1 for a power-law profile. 



5.4 Counter-rotating case 

In this Section, we investigate the counter-rotating configu- 
ration for a BH and its misaligned accretion disc, for initial 
values of O ut,o close to n. 

As shown by |Scheuer fc Feiler| |l996[ ) and |Martin et al.| 



the outer regions of the accretion disc, if this disc is regularly 
and coherently fed. Due to the Bardeen-Petterson effect, we 
expect the innermost part of the disc (approximately within 
7? warp ) to orbit in a plane which is perpendicular to Jbh, 
with orbital angular momentum density L counter-aligned 
with respect to the BH spin. In this BH-disc configuration, 
one of the major changes is in the radius of the innermost 
stable orbit, which increases due to the asymmetry seeded 
in the geodetic motion of particles in Kerr metrics. As a 
consequence, the energy and the orbital angular momen- 
tum of particles at Riso increase, while the BH radiative 



efficiency decreases (see, i.e., Wilkins 1972 



|1972[ ). The Bardeen-Petterson timescale (18 1 and the warp 



Bardeen et al. 



radius |19| have the same values as in the co-rotating case. 
Since tap <S t a i, the adiabatic approximation holds again, 
but the small deformation approximatiorj^jhas a limited va- 
lidity, requiring 9 ou t ~ n. In order to remain consistent with 
the approximation scheme, we trace the alignment process, 
from 7T to (it — tt/6), only. 

In the couter-rotating case and small deformation ap- 
proximation (i.e. # out ~ 7r), the disc profile can be solved 
analytically. We choose a reference fram e Q" x"y" z" where 
Jbh = (0, 0, —1) and we solved equation ( 14 1 for 1 in it. For 
constant viscosities, the function W" (R/ 'i? W arp) = W + ily 



(20071, on the timescale t a l the BH spin again aligns with 



2 The functions ip and x, defined in equations ( |2_'| and \30\ , arc 
invariant under the transformation 9 ou t — » (tt — 8 ou t). 
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Figure 12. Color coded map of the aligment time (defined as the time necessary for the relative inclination angle to go from O ut,o = t/6 
to flout = tt/600) for a BH of M BHj0 = 1O 6 M , as a function of the accretion rate expressed through the Eddinfton factor /Edd an d of 
the initial BH spin parameter arj. The colour scale represents t a i in years. Top (bottom) panels refer to the constant (power-law) viscosity 
profile. Left (right) panels refer a = 0.09 (a = 0.18). 



describing the warp is 
W^ cnt = C exp ( -V2 + i 



while for a power-law viscosity profile 

i 

R 



= D 



x K 1 



/2(l+/3) 



V2(l + i) 
(1 + /3) \R 



R 



(55) 



(56) 



We then apply the adiabatic approximation to study the 
coupled evolutions of the system BH-disc. The jointed evo- 



lutions of 8 ou t, Mbh and a are presented in Figure [15] and in 
Table [3] for different accretion rates and viscosity profiles. 
The shorter timescales in the counter-rotating configuration 
stem from the the dependence of the alignment timescale on 
the spin modulus, i a i oc a 5 ^ 7 . Counter-rotating matter car- 
ries larger and opposite angular momentum, reducing the 
spin modulus and the alignment timescale in the process. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we followed the joint evolution of the mass 
Mbh and spin Jbh of a BH inside a geometrically thin, ex- 
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panels refer to the constant (power-law) viscosity profile. Left (right) panels refer a = 0.09 (a = 0.18). 



tended accretion disc. The BH spin is initially misaligned 
with the angular momentum of the disc in its outer re- 
gions. On the short Bardeen-Petterson timescale, the disc 
responds to the Lense-Thirring precession, imposed by the 
BH spin, and propagates a warp that is maximum around 
fl walp ; within this radius matter orbits around the BH in a 
plane which is perpendicular to the BH spin. According to 
angular momentum conservation, the warped disc interacts 
with the BH spin and, on the longer alignment timescale, 
the BH aligns its spin to Jdisc.out- In its outer regions, the 
disc is assumed to be fed by matter that flows along a plane 
that keeps its coherence in direction Jdisc.out for a sufficiently 
long timescale to allow for the gravitomagnetic interaction 
to complete BH-disc alignment. While doing so the BH is 



accreting matter and angular momentum from the inner por- 
tion of the disc, which is aligned or anti-aligned to the BH 
spin. Given the mismatch between the timescale for warp 
propagation and the alignment time ( Scheuer fc Feiler|1996 



|Natarajan fc P ringlc 1998]) , we devised a method that en- 
abled us to follow, in the small-deformation approximation, 
the co-evolution of the BH mass and spin in a self-consistent 
manner, carrying out a large survey of the parameter space 
and a critical review of the used approximations. 

It is found that, considering an initial small relative in- 



clination angle (0 o 



7r/6, small deformation approxima- 



tion) , matter in the inner part of the accretion disc has or- 
bital angular momentum density parallel to Jbh- The gravit- 
omagnetic interaction of the BH with this warped accretion 
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Figure 14. Color coded map of the relative increase of spin parameter during the alignment time (defined as the time necessary for the 
relative inclination angle to go from 9 ou t : o = t/6 to 9 ou t = 7r/600) for a BH of Mbh o = 10 6 Mq, as a function of the accretion rate 
expressed through the Eddinfton factor /Edd an d °f the initial BH spin parameter ao- The colour scale represents Aa/ao- Top (bottom) 
panels refer to the constant (power-law) viscosity profile. Left (right) panels refer a = 0.09 (a = 0.18). 



disc and their coupled evolution bring the BH into align- 
ment with the outer regions of the disc, i.e. O ut(t) — > 0. 
The timescale t a i of equation ( 43 1 gives a good estimate of 



the BH-disc alignment time for an e-folding reduction of the 
angle of misalignment, in very good agreement with numer- 
ical results. For a maximally rotating Kerr BH accreting 
at the Eddington rate, t a i ~ 10 5 ~ 6 yr, depending on the 
viscosity parameter a and on the viscosity profile model, 
in agreement with early findings by |Natarajan fc Pringle| 
(19981. On the other hand, environments where the accre- 



tion rate is extremely low imply longer alignment timescales, 
as t a i oc M ~ 32 / 35 . I n the explored BH mass range, the align- 
ment time displays a weak dependence on Mbh.o: fixed all 
the other parameters, alignment of a 10 7 M s BH occurs, on 



average, at the same pace of a 10 Mq BH. The BH mass 
and spin modulus increase during alignment, but their frac- 
tional increases are modest. After surveying a wide param- 
eters space, we find that 0.1% < AM B h/M B h,o < 3% while 
the spin parameter increases by 0.5% < Aa/ao < 20%. 

Starting with an almost anti-parallel BH-disc configu- 
ration (0 O ut,o ~ i"), the orbital angular momentum density 
of the inner part of the disc is initially counter-aligned with 
respect to the BH spin. Nevertheless, the BH still tends to 
reduce the degree of misalignment (i.e. 8 out (t) decreases), be- 
cause of the nature of the gravitomagnetic interaction (see 
also |Scheuer fc Feiler||l996| |Martin et al.||2007[ ). The accre- 
tion of matter with opposite angular momentum at Riso 
decreases Jbh and a with higher rates, compared with their 
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Figure 15. Coupled evolution of the relative inclination angle (n — 6 ou t), BH mass M^n and spin parameter a for a counter-rotating 
disc. Solid (dashed) lines refer to constant (power-law) viscosity profile. Black lines to /Edd = 1, red lines to /Edd = 0.1, blue lines 
to /Edd = 0.01. Dotted horizontal lines which appear in top panels represent angles 0out/#out,o = 10, 10 2 ,10 3 . Initial configuration: 
M B h,o = 10 6 M o , a = 0.5, / Ed d = 0.1 and a = 0.09, with tt - ou t,o = tt/600. 



growths in the specular co-rotating case. Since t a i oc a 5 / 7 , 
the alignment process is then more efficient, i.e. the angle 
reduction speed is higher. Comparing decreases of the rel- 
ative inclination angle # out symmetric with respect to 7r/2, 
we find that the fractional decrease of the spin parameter in 
the counter rotating case is, in modulus, higher than in the 
specular co-rotating case while the mass relative increase is 
slightly lower. BH spin flip, due to 8 out reduction below 7r/2, 
will occur in this extended disc when a will reach its mini- 
mum value. At that time the jet of relativistic particles (if 
present) will cross the warped disc, likely affecting the sub- 
sequent BH-disc evolution and the BH feeding. This process 
deserves a separate investigation. 

It is still poorly known whether a spinning BH in an ac- 
tive galactic nucleus is fed through a disc that maintains its 
angular momentum direction stable over a Salpeter timescale 
ts- Two opposite, still plausible scenarios, have been pro- 



that the stability of jets in radio-loud AGNs requires a long- 
lived phase of stable accretion capable to maintain spatial 
coherence, i.e. a fixed direct ion of Jdisc.out, for a time a s 
long as 10 s yr. By contrast, |King fc Pringle| (|2006| |2007[); 



King et al. 



( 2008 1 speculated recently that AGN activity, 
is chaotic in nature 
is occurring through 



triggered by gas-rich major mergers, 
even within a single merger event, i.e. 
a sequence of uncorrelated short-lived accretion episodes. 
In their picture the corresponding discs, truncated by their 
own-self gravity, continuously change their inclination and 
feed the BH on their consumption timescale. Under these cir- 
cumstances the BH spin modulus is seen to either increase 
or decrease at random clustering around small average val- 
ues a ~ 0.1 — 0.3. This model would simultaneously explain 
the relatively low radiative efficiency of the quasar popu- 



lation as inferred from the background light (e.g., Merloni 
2004 Merloni & Heinz 2008), and the possibility of growing 



posed and discussed. Natarajan & Pringle ( 1998) speculated 
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Table 3. Summary of our parameters and results for the counter-rotating case; we consider viscosity coefficient a = 0.09 and initial 
inclination angle 9 ont ,o = 7r (l — 1/6000), both for constant (C) and power-law (PL) viscosity profiles (VP). The initial BH has A/bh,o = 
1O 6 M0 and aa = 0.5. Accretion rate /Edd varies over three orders of magnitude and we record times needed to (n — O ut, 0) of a factor 
10, 100 or 1000; we also report mass and spin relative variations. 



BH as massive as 10 Mq from small BH seeds already at 
redshift z ~ 6 ( |King fc Pringle|2006 1. 

Isolated discs, truncated by their own self-gravity, carry 
a well defined disc angular momentum and are accreted by 
the BH on a finite timescale. Starting with a misaligned 
BH-disc configuration, the BH spin changes direction sig- 
nificantly only if (i) the alignment time is shorter than the 
disc consumption time, t a i < tdi sc ; and (ii) the magnitude of 
the disc angular momentum is comparable to the BH spin 
magnitude, i.e. Jdisc > Jbh- The first condition is verified 
for the whole parameter range explored in this paper. The 
estimate Jdisc > Jbh depends instead sensitively upon i? ut- 
Equation[l2]establishes that isolated discs around large BHs 
truncate at -R ou t such that Jdisc < Jbh- Condition (ii) is sat- 
isfied for BH masses < 3 x 10 7 Mq. We note here that since 
we model our discs using the Shakura-Sunyaev solution for 
Kramer's opacity, we cannot rigourously estimate R ou t, and 
Jdisc, for BHs with mass < 10 5 — 1O 6 M0. An extension of disc 
solutions to different, self-consistent opacities is non trivial 
I Hure et al.|1994a|b I and we postpone a detailed analysis to 
future work. 

BHs with masses < 3 x 10 7 align efficiently in discs 
truncated by their own self-gravity, implying alignment also 
in the case of stochastically fed AGNs^l where not only a 
fluctuates with time, but also the direction of the BH spin 



3 Here we are not considering accretion events involving a disc 
with a very small amount of mass, i.e. below ~ (H / R)Mbb\- This 
light accretion disc has an outer radius much smaller than (JsJ 
and thus carries an angular momentum Jdisc *C Jbh- As a con- 
sequence, the disc has a very short consumption timescale, and 
the alignment process is active for a very short period of time. 



continually changes due to the rapidity of the alignment 
process. By contrast, rapidly spinning (a ~ 1) heavier BHs 
with Mbh ^ 1O 8 M0 have truncated discs that carry little 
angular momentum compared with Jbh- In this case align- 
ment is uneffective and the orientation of the BH spin is not 
influenced significantly by the surrounding short-lived disc. 

In light of these findings, the vector Jbh appears to 
carry precious information on the orientation of the plane 
through which the BH has been fed, and on whether ac- 
cretion has been long-lived and coherent or short-lived and 
random. 

The method developed in the paper is sufficiently ver- 
satile that it will be implemented in numerical simulations 
describing the process of pairing of dual BHs in circum- 
nuclear discs during their on-fly accretion (Dotti et al., in 
preparation) to improve upon the speculation ( |Bogdanovic| 
|et al.||2007| ) that, in gas-rich galaxy mergers, binary BHs 
have time to align their spin orthogonally to their orbital 
plane, as discussed in Escala et al. (120051) ;|Dotti et al. 1(2006); 



|Mayer et al] ( [2007] ; [Dotti et al.| ( |2007H2009[ ); |Colpi fc Dotti 
(|2009h. The spin-orbit configuration is relevant to study the 
impact of BH recoils, that occur after two BHs have coa- 
lesced (see, e.g., Pretorius||2007 1 . Detection of gravitational 
waves, emetted by coalescing BHs, with the Laser Interfer- 



ometer Space Antenna (LISA) (Bender et al. 1994 Hils & 
|Bender||l995[ ) will be able to constrain the moduli and the 



directions of the coalescing BHs spins ( |Vecchio|[2004l |Lang 
fc Hughes|2006 i. 



Therefore the alignment of Jbh around J to t (close to Jbh) is 
expected to be unimportant. 



20 A. Perego et al. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF 
THE GRAVITOMAGNETIC TORQUE 

Constant viscosity: For the constant viscosity model, the 



disc profile is described by ( 29 \ and the equation ( 53 1 be- 
comes 



S(Jbh,x + iJBH,y) 
2V2A 



(i-0- 



x exp 



V GMbH : — tep(^warp) 



A r 2 R 5/4 



(Al) 



-72(1 



R 



Rv 



-1/2 



Power-law viscosity: For the power-law case Q with ex- 
ponent /3, the disc profile is given by Wpl, defined as (281 



(2007) 



In this case, equation ( 53 1 was integrated by Martin et al 



gin 



x t B p(R warp 



8GM JbhV GMbh 

4/3 + 3 

JLwarp 



2 °ut 20 + 1 

Z 2(T+P)K i (z) dz 

2(1 + 0) 



(A2) 



where z is a new complex variable, defined as 



^2 



R 



Rv. 



(A3) 
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Assuming that 

Zout 2/3 + 1 

zW+mK i (z) dz « 

in 2(1 + 0) 

(l-i)oo 2 g+l i 

z 2(1+0) K_i_(z) dz = 2 str+W (A4) 

2(1 + 3) 

1 + 2/? 



ii 



x r 



2(1 + P) 



we can rewrite the infinitesimal gravitomagnetic spin varia- 
tion as 



8(JbH,x + iJ'BH,j/)gm _ ■- 4(i + fl) ^BP(-Rwarp) 

where 

1 _ 4GMVGMBH 7 R / V2 



2/3 + 1 
' 2(1 + 0) 



2(1 + 



(A5) 



(A6) 



Martin et al. (20071 estimate the alignment timescale as 



Ui,u = 7^ r - (A7) 



(4(14-/8)) 



